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Abstract 

The maximum complexity momentum distribution for an isolated monodimensional ideal gas 
out of equilibrium is derived analytically. In a first approximation, it consists of a double non- 
overlapping Gaussian distribution. In good agreement with this result, the numerical simulations 
of a particular isolated monodimensional gas, which is abruptly pushed far from equilibrium, shows 
the maximum complexity distribution in the decay of the system toward equilibrium. 
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1. INTRODUCTION 



Boltzmann-Gibbs (BG) statistics works perfectly for classical systems in equilibrium 
under the action of short-range forces. But most systems in nature are out of equilibrium and 
there is no a priori reason why a particular phenomenon should behave according to a specific 
kind of statistics. If the system finally decays toward equilibrium, then the asymptotic long- 
time limit should be that of BG statistics. 

Different kinds of statistics have been proposed to model different nonequilibrium sit- 
uations. For instance, non- extensive thermo statistics is based on maximizing the Tsallis 
entropy |j| under different assumptions for calculating the expectation value of the energy. 
Power law distributions are obtained by fixing the total energy of the system in all the cases 
analyzed in Ref. ^. Thus, in this scheme, the exponential distribution of BG statistics 
turns out to be a singularity that is recovered in the limit g ^ 1, where q is called the index 
of non-extensitivity. Although this type of statistics might seem a mathematical artifact 
without applications, several types of generahzed stochastic dynamics have been recently 
constructed for which Tsallis statistics can be proved rigorously P]. Also, it has been found 
useful in explaining many other physical phenomena jj, O]. Another formalism to study 
out of equilibrium situations is super statistics. Originally proposed by Beck and Cohen 
0], it deals with nonequilibrium systems with a long-term stationary state that possess a 
spatio-temporally fiuctuating intensive quantity. After averaging over the fiuctuations one 
can obtain an infinite set of general statistics called superstatistics, which constitute a su- 
perposition of BG distributions. Tsallis statistics is a special case of such superstatistics, 
and, in particular, BG statistics is also recovered when g — * 1, where q is now a dynam- 
ical parameter with a certain physical interpretation. In general, complex nonequilibrium 
problems may require different types of superstatistics :6]. 

Although all these statistical techniques for modeling out of equilibrium situations are 
becoming a well established theory, as far as we know, there are no general laws telling us in 
what manner a system should relax towards equilibrium. The second law of thermodynamics 
claims that the average entropy or disorder must increase when an isolated system tends to 
equilibrium but no more insight is obtained from this postulate. In fact, this law in no way 
forbids local complexity from arising An inspiring example is life, which can continue 
to exist and replicate in an isolated system as long as internal resources last. It could then 
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be postulated that in an isolated system, besides an increase in entropy, the system will try 
to stay close to the maximum complexity state. This behavior was found in Ref. ^ for a 
particular system, the "tetrahedral" gas, when complexity is defined as in Ref. 0] (referred 
to in the literature as the LMC complexity). Furthermore, it was established that this 
isolated system relaxes towards equilibrium by approaching the maximum complexity path. 
This path is an attractive trajectory in the distribution space connecting all the maximum 
complexity distributions (see Ref. j8| for details). 

In this paper, we perform the study of an isolated mono dimensional ideal gas that is 
initially in equilibrium, receives a strong perturbation, and finally freely decays towards 
equilibrium. The one-particle momentum distribution is computationally calculated for 
each time during the relaxation process. It is found that this distribution coincides with 
the maximum LMC-complexity one-particle momentum distribution. Hence, the maximum 
complexity path in the space of one particle momentum distributions also seems to explain 
the statistical evolution of this system when it approaches equilibrium. 



2. NUMERICAL SIMULATIONS OF A MONODIMENSIONAL GAS 

In order to have a graphical picture of the processes involved, we start by describing the 
numerical simulations. The gas is initially, at time t = 0, in equilibrium. Its one particle 
momentum distribution is described by a Gaussian or Maxwell-Boltzmann function. At this 
point, two new extremely energetic particles are introduced into the gas, forcing the gas into 
a far from equilibrium state. The system is kept isolated from then on. It eventually relaxes 
again toward equilibrium showing asymptotically another Gaussian distribution. Most of 
the time, during this out of equilibrium process, the momentum distribution function is 
described to a first approximation by two Gaussian distributions. This double Gaussian 
distribution coincides with the analytically derived maximum complexity distribution, which 
will be derived later. 

In more detail, and using arbitrary units from now on, the gas consisted of 10 000 pointlike 
particles colliding with each other elastically. The particles were positioned with alternating 
masses of 1 and 2 on regular intervals on a linear space 10 000 units long. The system has no 
boundaries, i.e., the last particle in this linear space was allowed to collide with the first one, 
in a way similar to a set of rods on a circular ring. Two distinct masses in the system were 
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FIG. 1: Results of the numerical simulations of an isolated monodimensional ideal gas when it 
relaxes towards equilibrium (lower right) from an initial condition in which two very energetic 
particles are introduced into the gas in equilibrium (upper left). Histograms of the one-particle 
momentum distribution with a log vertical axis at various times are shown as solid lines. A double 
Gaussian fit, the maximum complexity distribution, is shown as dashed lines. 



used because a monodimensional gas can thermalize only if its constituent particles have at 
least two different masses. Initially 9998 particles were given initial conditions following a 
Gaussian distribution with mean zero velocity and a mean energy of 1/2, giving a total mean 
energy for the system of nearly 5000. These particles where then allowed to undergo 20 mil- 
lion collisions in order for the system to reach the initial state of equilibrium, i.e., a Gaussian 
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distribution. After that, at time t = two extremely energetic particles of mass 1 and 2 
are introduced at two neighboring points, such that the total system has zero momentum 
and an energy of 150 000. The system then undergoes another 20 million collisions to reach 
again the equilibrium Gaussian distribution after a total elapsed time of At = 1156.05. We 
record the time evolution of the one- dimensional momentum distribution in Fig. ^ where 
the square of the generalized particle momentum is given by the variable pf = pf/rrii, with 
Pi and rrii the momentum and mass of particle i. The theoretical maximum complexity dis- 
tribution (derived below) for this system is approximated by two non- overlapping Gaussian 
distributions and is also fitted as dashed lines in Fig. ^ Let us remark that the system 
stays in this double Gaussian, the maximum complexity distribution, during a large part 
of its out of equilibrium phase. The two clearly visible slopes of Fig. ^ are related with 
the two different mean energies associated with both Gaussian distributions. As the system 
approaches equilibrium, both Gaussian distributions merge into one. 

In Fig. 121 the momentum distribution is shown at a particular moment out of equilib- 
rium (solid line). In this case a direct histogram of the momentum distribution is shown. 
The double gaussian or maximum complexity distribution is clearly seen. In this Figure, 
the complete gaussian distributions are shown for clarity (dashed lines), but the fit to the 
numerical simulated data has been done with non-overlaping Gaussian functions. 



3. MAXIMUM COMPLEXITY DISTRIBUTION 



We now sketch the analytical derivation of the microcanonical maximum complexity 
distribution for a system with a huge number of accessible states. In Ref. |8|, the maximum 
complexity distribution for an isolated system with a discrete number of accessible states was 
derived. This type of distribution was found to be important in the path towards equilibrium 
for a particular isolated gas, the tetrahedral gas. When this system is out equilibrium it 
decays to equilibrium by approaching the trajectory formed by all the maximum complexity 
distributions and called the maximum complexity path. Finding these extremal distributions 
requires solving a variational problem. The complexity C (see Ref. Q), is defined as 

C = D-H, (1) 
where the disequilibium, D, is defined as the distance to the microcanonical equilibrium 



5 



•f collisions: 2. 



Time: 202.51 



m 

1j 



o 



o 
o 
^ 

o 
o 

CO 

o 
o 

o 
o 

^ — I 

o 




20 







20 



P 



n 



FIG. 2: Numerical simulations of an isolated monodimensional ideal gas as it relaxes to equilibirum 
at a particular time, 202.5 after 2 million collisions. The histogram of the one particle momentum 
distribution is shown as a solid line. The two fitted gaussian distributions, the maximum complexity 
distribution, are shown as dashed lines. 

probability distribution, the equiprobability, and H is the normalized entropy, 
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D = Y,{U - l/Nf and if = -(1/ In N) J] /, In /, 



(2) 



i=l 



where is the number of accessible states and /j, with i = 1,2, . . . , N, are the proba- 
bilities of permanence that the system presents for the different discrete accessible states i. 
Thus, the micro canonical maximum complexity distribution can be derived by finding the 
maximum disequilibrium for a given entropy using Lagrange multipliers . The results are 
shown in Table |l] Note that, for a given entropy, maximizing the disequilibrium is equivalent 
to minimizing the Tsallis entropy with parameter q = 2. Note also that for an isolated sys- 
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tern the entropy variable is equivalent to a streched time scale, due to its monotonic increase 
with time given by the second law of thermodynamics. The results of Table HI graphically 
represented in Fig. El show that the maximum complexity distribution can be split into 
two components. One of them consists of a background equiprobability distribution for all 
accessible states, which will be denoted as the "people distribution". The other one, with 
the remaining probability, comprises the particular state with the highest probability and 
is called the "king distribution" . The final maximum complexity distribution is the sum of 
the people and the king distributions. When the system reaches the equilibrium, the king 
and people distributions merge leaving only the equiprobability distribution. 

4. MAXIMUM COMPLEXITY DISTRIBUTION OF AN ISOLATED MONODI- 
MENSIONAL IDEAL GAS 

The derivation of a maximum complexity distribution will be based on the symmetry 
of the momentum phase space and on particular initial condition considerations. As an 
extension of the former discrete case given in Table lU the distribution we are looking for 
will have two components. The people distribution component will be the equiprobability 
distribution and the king distribution will be made by choosing a particular dominant state 
according to plausible arguments derived from the initial state of the system. 

First, let us obtain by symmetry arguments how the functional dependence of the equi- 
librium distribution, or equivalently the equiprobabilty one, looks when only one variable 
is maintained and all the rest are integrated. In an isolated ideal gas with n particles, all 
accessible states lie on the surface of a hypersphere in the p„-momentum phase space. If the 
energy of a single particle is Cj = then the total energy of the ensemble is E = Yl^=i 
The mean energy per particle e will he e = E/n. If the gas is in equilibrium, the micro- 
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FIG. 3: Micro canonical maximum complexity distribution derived in Ref. 1^ (thin solid line). It 
is the sum of the king (thick solid line) and the people, or equiprobability, distribution (dotted 
line). 

canonical distribution, h, is the equiprobability for all accessible states, i.e., all points in 
phase space lying on the hypersphere surface have the same weight in the distribution h. 
This distribution is given by the expression 

h{9n-l, On~2, ■ ■ ■ , 9i)d9n~ld9n-2 ■ ■ ■ d9i = 

r'^~^d9n-i sin 9n-id9n-2 sin 9^-1 sin 9n~2d9n-3 ■ ■ ■ 

sin 9n-i sin 9n-2 ... sin 92d9i , (3) 

where the original phase space variables, {pi,p2, . . . ,Pn), have been converted to the spherical 
coordinates, {r,9i,92, ■ ■ ■ ,9n-i), with = ne. To obtain the one-particle momentum 
distribution of this system in equilibrium h can be integrated over all coordinates except 
9n-i obtaining the function g, 



g{9n-i)d9n-i = C'r^-'siYi^-'9n-id9, 



n— 1 ; 



(4) 
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where C is the constant of integration. Converting this result back to the ]3„-momentum 
coordinate via the relation p„ = r cos^^„_i we have the final distribution /, 

n — Z 

f{Vn)dVn = C'r^^' (1 - Pl/r') ^ dp^. (5) 

Taking the limit for a large number, n, of particles the one-particle momentum distribu- 
tion in equilibrium is obtained, 

fEQ{Pn)dpn = Cexp{-pl/4:e)dpn, (6) 

where C = y^TJjAnej is the normalization constant. Substituting e = KT/2, K being the 
Boltzmann constant and T the temperature, the familiar one-particle ideal-gas momentum 
distribution is obtained, i.e., the Maxwell-Boltzmann distribution. It is remarkable that this 
result has been obtained in the microcanonical ensemble although this distribution is usually 
presented as a typical derivation from the canonical formalism. The people distribution, 
/p, being an equiprobability distribution, its formal expression will be the same as the 
equilibrium one, 

fp{Pn) = Cp exp(-p2/4ep). (7) 

Fig. m shows the people distribution in a dimensional space with just three particles as a 
dashed mesh surface. 

Second, let us discuss how to find a king distribution for our system. In the scenario of 
injecting two extremely energetic particles into the gas in equilibrium, the states with the 
highest velocities will be populated with a certain number of particles in phase space. Using 
symmetry considerations, all high-velocity components in phase space for each one of those 
particles should be equivalent in this distribution. It then seems plausible to assume that 
the high-probability state of the king distribution can be formed by all the spherical caps in 
phase space centered on the maximum possible velocities for each particle. Fig. E] illustrates 
this concept in a small dimensional space with just three particles. The king distribution in 
this phase space is shown as a mesh of solid lines: the probability of a state lying inside one 
of these caps is uniform and zero if it lies outside. Each of these caps is identified by the 
angle a shown in the figure. Integrating for all momentum dimensions except one, the one 
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FIG. 4: Non-zero constant probability surfaces of the people (dashed mesh) and king (solid mesh) 
distributions on the space {3D sphere) of accessible states of a three particle monodimensional 
ideal gas. 

particle momentum distribution Jk is obtained, 



$(p„) ifp„<rsinQ;, 
fK{Pn) = < if r sin a < Pn < r cos a, (8) 

Ck e*^"^"/"^*^^^ if r cos a < Pn, 
Since the spherical caps that are over the p„ = point only cover a part of the hypersphere 
surface given by the set of points (pi,P2, ■ ■ ■ ,Pn-i,Pn — 0), the function $(p„) will satisfy, 
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< CKexp{-pl/AeK). 

The final exact expression for tlie maximum complexity non-equilibrium distribution will 
be the sum of both distributions, 



fMciPn) = fKiPn) + fp{Pn)- (9) 

We can simplify this equation by noting that, for the Pn < r sina cases, $(pn) ^ fp{Pn), 
and that, for the rcosa < p„ cases, fpipn) ^ fK{Pn)- 

As a first approximation, the maximum complexity non-equilibrium distribution is ob- 
tained: 

\Cpexp{-pl/4ep) if p„ < po, 
fMciPn) = < (10) 
[^C^exp(-p2/4e^) if < p^^ 

where the momentum pQ = r cos a is defined. Hence the final maximum complexity 
distribution is a two-component non-overlapping Gaussian function characterized by the 
parameters ep and ex, which vary with time when the system decays towards equilibrium, 
as seen in Fig. ^ 



5. CONCLUSION 



In conclusion, a maximum complexity one-particle momentum distribution has been de- 
rived for an isolated monodimensional ideal gas far from equilibrium. It is based on maximiz- 
ing the disequilibrium, or equivalently minimizing the Tsallis entropy with parameter q = 2, 
for a given entropy, or equivalently time, in an isolated system. In a first approximation, 
the maximum complexity distribution is a double non-overlapping Gaussian distribution. 
Numerical simulations of a particular isolated monodimensional gas show, in clear agree- 
ment with our analytical result, a double Gaussian distribution when it decays towards 
equilibrium. 
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